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^ ■ ABSTRACT 

We investigate the amount of scatter in the Tully-Fisher relation 
(TFR) when using optical long-slit Ha rotation curves to determine the 
velocity widths of spiral galaxies. We study a sample of 25 galaxies in 
the Coma region of the sky which were shown in Bernstein et al. (1994) 
to exhibit an extraordinarily low scatter of 0.10 mag RMS in the / mag- 
nitude vs 21-cm width TFR. Using the same / magnitudes with new 
widths derived from high-quality Ha rotation curves, we measure an 
RMS scatter of 0.14 mag in the TFR. This suggests that measurement 



1 Based on observations obtained with the Multiple Mirror Telescope, a joint facility of the Smithsonian Institution 
and the University of Arizona. 
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errors and "astrophysical errors" (such as non-circular gas motion) on 
the Ha velocity widths are below 6%, and optical widths are nearly as 
good for TFR studies as 21-cm widths. The scatter and form of the 
TFR are found to be robust under choice of velocity width-extraction 
algorithm, as long as the central portions of the optical rotation curve 
are ignored and low-S/N points are not weighted too heavily. In this 
small sample there is no evidence that rotation curve shapes vary sys- 
tematically with rotation velocity, nor that rotation curve shape can be 
used to reduce the scatter in the TFR. 

Subject headings: cosmology: distance scale; galaxies: spiral; galaxies: 
fundamental parameters; galaxies: distances and redshifts 

Accepted for publication in the Astronomical Journal 
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1. Introduction 

Peculiar velocity surveys of thousands of 
spiral galaxies are now being conducted us- 
ing the Tully-Fisher relation (TFR, Tully & 
Fisher 1977) between the speed of rotation of 
a spiral and its absolute magnitude. The ma- 
jority of these surveys have been conducted 
using 21-cm Hi rotation widths, but many 
now also make use of long-slit optical Ha ro- 
tation curves to determine the rotation speed 
(e.g. Matthewson, Ford, & Buchhorn 1992, 
Giovanelli et al. 1997ab, Willick et al. 1996). 

In this paper we ask how the scatter in 
the TFR using optical data compares to the 
21-cm TFR scatter. This comparison is fa- 
cilitated by the use of a sample of 25 galax- 
ies in the Coma Supercluster region which we 
have previously shown to have extraordinar- 
ily low scatter in the I magnitude vs 21-cm 
width TFR — only 0.10 mag RMS (Bernstein 
et al. 1994, hereafter referred to as Paper I). 
We have obtained high-quality Ha rotation 
curves for all these galaxies, have extracted 
velocity widths, and determine the TFR scat- 
ter using optical widths. Since the photom- 
etry and 21-cm data for these galaxies (ei- 
ther by chance or for some not-yet-understood 
physical reason) give a nearly perfect 21-cm 
TFR, we are able to perform a very sensitive 
test for additional scatter induced in the TFR 
by a switch from 21 cm to Ha linewidths. The 
scatter in the optical TFR will present an up- 
per limit on the TFR scatter associated with 
the optical widths, whether the cause be mea- 
surement error or some failure of the rotation 
curves to reflect the "true" rotation width. 

2. Expectations and Previous Results 

A priori, 21-cm widths are expected to be 
superior to optical widths in the context of 



TFR studies. The detectable 21-cm emission 
extends further out into the disk of the spiral 
than the Ha emission, and the gas motions 
are likely more regular and circular in the 
outer disk than in the central sections. The 
21-cm velocity should then better reflect the 
presumed underlying physical variable of the 
TFR, galaxy mass, and thus produce lower 
TFR scatter. Secondly, all Hi atoms in the 
galaxy contribute equally to the 21-cm pro- 
file (dubbed "atomic democracy" by Schom- 
mer et al. 1993) so the data are not biased 
or restricted to streets of high ultraviolet ra- 
diation. Furthermore, 21-cm signals are not 
attenuated by dust — though extinction of Ha 
photons is relatively small at 2-3 scale lengths 
into the disk where the most relevant dynam- 
ical information lies. Finally, it is easy to 
achieve < 10kms~ resolution in 21-cm ob- 
servations. 

Long-slit rotation curves have the advan- 
tage over 21-cm profiles of being a two-dimen- 
sional data set (spatial and velocity axes). 
The latter have no spatial information, so 
it is not possible to distinguish thermal or 
non-circular motions from true rotation. On 
the long-slit rotation curve, random motions 
appear as noise atop the true rotation and 
thus can be removed. In choosing an algo- 
rithm to extract a single rotation width from 
the rotation curve, we may also restore some 
"atomic democracy" if we do not place too 
much weight on Hn regions with high flux, 
even though such regions might have very 
small formal errors on the recession velocity. 
Furthermore the morphology of the rotation 
curve could provide extra information to help 
reduce TFR scatter. 

It is also possible to obtain 3-dimensional 
data (two spatial axes and one velocity axis) 
using Fabry-Perot images of the Ha emission, 
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as demonstrated by Schommer et al. (1993). 
In theory this extra information should help 
detect dynamical irregularities even better 
than long-slit data, and also allows deter- 
mination of the galaxy inclination and po- 
sition angle independent of the photometric 
data. By simulating long-slit observations 
of their Fabry-Perot datacubes, Schommer et 
al. demonstrate several potential problems for 
long-slit measurements. They conclude, how- 
ever, that for highly inclined galaxies differ- 
ences between Fabry-Perot and long-slit esti- 
mates of the rotation speed will be small. 

Long-slit Ha measurements are perhaps 
the easiest to make in practice. They cer- 
tainly require less telescope time than Fabry- 
Perot measurements, and can be extended 
to much more distant galaxies than 21-cm 
methods — especially in regions of the sky which 
are inaccessible to the Arecibo telescope. Ro- 
tation curves of spirals at 7000 km s -1 reces- 
sion velocity can be obtained in under an hour 
with 2-meter class optical telescopes, but be- 
come difficult to observe with radio telescopes 
other than Arecibo. Of course a practical dis- 
advantage to long-slit observations is that one 
must know the position angle of the galaxy 
dynamical major axis in advance (preferably 
to an accuracy of 2°-3°). In our study the 
photometric position angle is known to 1°- 
2° because we precede the Ha observations 
with very deep /-band imaging, and we re- 
move from our sample any galaxy in which 
the outer isophotes twist by more than a few 
degrees. It has been noted (Franx & de Zeeuw 
1992) that if spiral disks are not intrinsically 
round, then the photometric major axis will 
not in general coincide with the dynamical 
major axis, causing errors in long-slit observa- 
tions. Our optical TFR scatter may be taken 
as a test of this non-circularity. 



These practical considerations have led oth- 
ers to attempt TFR measurements with op- 
tical rotation curves. Courteau (1992) con- 
ducted a TFR peculiar velocity survey with 
Ha rotation widths, and investigated several 
algorithms for extracting a width W from a 
rotation curve. These are evaluated first by 
checking variations on repeat measurements, 
and also by comparison with 21-cm values. As 
one might expect, algorithms which involve 
choosing a maximum point on each arm of 
the rotation curve fare poorly because they 
will usually extract the width from the nois- 
iest bins. The lowest internal scatter in W 
(lOkms -1 ) results from essentially averaging 
all velocities over each arm of the rotation 
curve. Courteau chose, however, to add all 
the spectra along the slit, discarding all spa- 
tial information, and collapsing all the Ha 
flux into a single velocity profile. He then 
derived W by means similar to those used for 
21-cm profiles. While his spatially-averaged 
optical widths gave higher internal errors (12- 
13kms _1 ) than W estimates which use the 
spatial information in the rotation curve, the 
offset from 21-cm H^'s is smaller, which made 
them better suited to his purposes. 

Other investigators have used Ha rotation 
curves to obtain W but do not explicitly com- 
pare various algorithms. Mathewson, Ford, 
& Buchhorn (1992) subtract the minimum 
from the maximum of the rotation curve to 
obtain W, and measure an internal scatter 
of 10 kms -1 . As discussed below, it is likely 
that they have enforced some unspecified min- 
imum S/N on the rotation curve points. Gio- 
vanelli et al. (1997ab) reinterpret some of 
the Matthewson, Ford, and Buchhorn rota- 
tion curves by defining W to be the width at 
R opt (the 83%-light radius) from the galaxy 
center. This requires extrapolation for some 
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rotation curves which do not reach this deep 
into the disk (3.2 scale lengths). 

Vogt (1994) extracted widths from rota- 
tion curves using a variant of the velocity- 
profile method: she measured the mean ve- 
locity in each spatial bin along the rotation 
curve, ranked these velocities in order, and 
defined the rotation width to be the difference 
between the 10th and 90th percentile points of 
the velocity distribution. This should encour- 
age more atomic democracy than the Courteau 
algorithm, because it weights all parts of the 
rotation curve equally rather than by flux. 

3. Measurements of Rotation Curves 
3.1. The Sample 

The Coma Supercluster galaxy sample and 
the collection of photometric data are de- 
scribed in Paper I. Briefly, we select those 
galaxies within a few degrees of the Coma 
cluster core, with recession velocities in the 
range 5000-8000 km s" 1 , for which 21 -cm pro- 
files could be found in the literature. Deep 
/-band surface photometry was obtained for 
each galaxy. We discarded galaxies for which 
the ellipticity varied by more than 0.03 across 
the outer isophotes or for which the position 
angle twisted by more than 3° — e.g., those 
with tidal tails, those for which the isophotal 
shape was dominated by the spiral arms, or 
those with morphological peculiarities. For 
these discarded galaxies (~15% of the full 
sample), we have little confidence that the 
isophotal shape accurately reflects the disk in- 
clination, making it difficult to derive the ro- 
tation speed from line-of-sight velocities along 
the photometric major axis. We also removed 
from the sample those galaxies (~6% of the 
full sample) for which the 21-cm profile did 
not have sufficiently steep sides to obtain a 



reliable width. 

In Paper I we demonstrate that few (if any) 
of the galaxies in our sample are likely to 
be members of the Coma cluster itself, and 
we henceforth assume that their distances are 
proportional to their redshifts. 

3.2. Observations 

Each of the galaxies was observed in March 
1993 or June 1993 using the Red Channel 
spectrograph on the Multiple Mirror Tele- 
scope (effective aperture of 4.5m). A 1200 
lines mm -1 grating blazed at 575 nm gives a 
resolution of ~0.21 nm per pixel near the Ha 
line. In addition, the thinned 800 x 1200 CCD 
was binned to 0'/6 pixels in the spatial direc- 
tion. The slit dimensions were l'/25 x 180", 
with the width chosen to admit as much light 
as possible without significantly degrading 
the spectral resolution. Most of the galaxies 
were observed for a single 1200 s exposure, 
with some of the fainter galaxies observed for 
two 1200 s exposures. This is significantly 
deeper than most existing TFR rotation curve 
studies — we wish to investigate the "intrin- 
sic" errors in the TFR using rotation curves 
rather than have the scatter be dominated by 
errors from photon noise. Position angles for 
the galaxies were chosen to match the outer 
isophotes of the /-band images. Galaxy nu- 
clei were centered on the slit using the MMT 
guiding camera. 

The galaxy observations were interspersed 
with short wavelength calibration exposures 
of a He-Ne-Ar arc lamp, and quartz lamp flat- 
field exposures. Quartz lamp spectra were 
also obtained through a "decker" slit consist- 
ing of a row of small holes. Long-slit spectra 
of the twilight sky were also obtained. 
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3.3. Extraction of Rotation Curves 

Each galaxy exposure was flat-fielded using 
the quartz lamp exposure adjacent in the ob- 
serving sequence, and a wavelength solution 
for the entire CCD was derived from the adja- 
cent He-Ne-Ar arc lamp exposure. Residuals 
to the wavelength fit are typically 0.005 nm. 
Twilight exposures were used to compensate 
for slight differences (dependent on the posi- 
tion of the slit) between lamp throughput and 
sky throughput caused by small differences 
between the pupil illumination by the lamp 
and by the sky. Finally, the decker slit quartz 
lamp exposures were used to map the distor- 
tion of the spatial coordinate on the chip. 
Thus for each exposure we produce a map 
from pixel coordinates (x, y) to (s, A) space 
(slit position vs wavelength). 

For a given pixel corresponding to (sj, 
we determined the sky background as follows: 
we determined the flux at Aj for each other 
row (s value) in the image. Note that since 
the spectrum is tilted relative to pixel coor- 
dinates, this may require interpolation. We 
then fit a linear function of s to the derived 
sky background intensity vector, omitting of 
course the central parts of the slit over which 
the galaxy emission can be seen on the CCD 
image. 

Cosmic rays were identified by eye, and the 
IRAF task IMEDIT was used to linearly in- 
terpolate over pixels affected by cosmic rays. 
Each pixel was then assigned to a bin accord- 
ing to its spatial coordinate s. Bin widths 
were 3 to 10 pixels (1//8-6") depending upon 
the S/N of the emission lines. Within each 
spatial bin, the intensities in each spectral 
pixel were fit to a model function of wave- 
length A which includes continuum radiation 
and gaussian profiles for the Ha and the two 
[N ii] lines (undersampling of the spectral line 



profile is accounted for). Each spectral line 
was allowed to have an independent intrin- 
sic width (which was then convolved with 
the instrumental width) and flux, but the 
three lines were constrained to have the same 
recession velocity. Thus for a linear back- 
ground, the spectral model has 2 parameters 
for continuum, 6 line strength/width parame- 
ters, and one velocity parameter. This model 
was fit to an ps 80-pixel (~17 nm) segment of 
the spectrum spanning the Ha and [N ii] lines, 
leaving pa 70 DOF for the fit. Chi-squared 
values were typically 70-100, except near the 
galaxy nucleus, where the continuum has high 
enough S/N to make the linear approximation 
inadequate. Note that we do not rebin or in- 
terpolate the spectra. 

The fitting process resulted in an estimate 
of recession velocity i>j and its uncertainty o~ Vi 
for each bin, with the bins centered at po- 
sitions Si along the slit. We only used bins 
in which the Ha emission was detected at 
S/N > 3.5, which leads to velocity uncertain- 
ties smaller than 12 km s _1 even in the noisiest 
parts of the rotation curves. Rotation curves 
for the 25 galaxies are shown in Figure 1. 

4. Velocity Width Algorithms 

We must extract a single width W from the 
rotation curves in order to produce a TFR. 
Desirable qualities for the W measurement al- 
gorithm are: 

1. The resultant W should be robust under 
changes in binning of the rotation curve 
data (or changes in spatial resolution). 

2. W should also vary little with the S/N 
level of the observations. 

3. W should have minimum variance for 
repeat observations at a given S/N level 
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Fig. la. — Rotation curves for the 25 galaxies in the sample. Bins with S/N below 3.5 are omitted. 
The x axis is in units of each galaxy's half-light radius ri/2, and the y axis is recession velocity 
in kms -1 . All galaxies are on a common scale with each vertical tick being 100 kms -1 . Dashed 
lines are the velocities of the two arms as determined by the Probable Min-Max method. 



of the spectrum. 

4. The algorithm should try to promote 
"atomic democracy" (see §2). 

5. The algorithm should not make a pri- 
ori assumptions about the shape of the 
rotation curve. 

6. Most importantly, it should minimize 
the scatter in the TFR. 



The first three of these criteria lead us to 
design algorithms that use as much of the ro- 
tation curve as possible in some form of aver- 
aging or fitting process — though criterion (5) 
leads us to avoid fitting functional forms to 
the curves. Choosing the maximum and min- 
imum velocity of the rotation curves (and sub- 
tracting to get W) is a very poor method, be- 
cause the final W is determined by the nois- 
iest bins. The derived W will have a strong 
tendency to increase as S/N worsens or spatial 
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Fig. lb. — Rotation curves for the 25 galaxies in the sample (continued). 



bin size is decreased. We describe here several 
width-determination algorithms which we ap- 
ply to the data in order to see which performs 
best under criterion (6), minimal TFR scat- 
ter. 

4.1. Min-Max Method 

This is the "straw man" method upon 
which we hope to improve. The width W 
is simply v max - v min , where v max and v min 
are the maximum and minimum of the i>j. 
This indeed leads to quite undesirable results, 
with a TFR scatter of 0.4 mag. The noise in 
the estimate of W is much abated if we con- 



sider only those V{ for which the Ha line has 
S/N > 3.5 and cr Vi < 12kms _1 ; we will hence- 
forth retain this restriction. Some form of 
this algorithm is used by Matthewson, Ford, 
k Buchhorn (1992), and Courteau (1992) 
also experiments with variants of this method. 
The "percentile" method of Vogt (1994)— 
choosing the width to be the difference be- 
tween the 10th and 90th percentile points of 
the velocity histogram — is similar to the min- 
max method, but with much reduced sensi- 
tivity to the noise level. Vogt also reduces 
sensitivity to noise by using only points with 
small uncertainties. 
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4.2. Weighted Mean Method: Entire 
Side 

The width W is here defined to be |i>i e ft — 
fright |, where the left and right side velocities 
are the weighted averages of the entire side 
of the rotation curve. The center is located 
at the maximum of the continuum flux. The 
velocity data points are weighted by <x~ 2 as 
usual; this might subvert "atomic democracy" 
by weighting too heavily those bright H n re- 
gions with high flux and thus low a Vi . We 
could encourage atomic democracy by replac- 
ing a v by some minimum value (~ 1 kms -1 ) 
if the flux is high enough to make the for- 
mal error smaller. For our sample, however, 
we find that this does not make a significant 
difference for any of the weighted mean algo- 
rithms considered. 

4.3. Weighted Mean Method: Outer 
Segments 

This algorithm is the same as the previ- 
ous algorithm, except that the inner part of 
the rotation curve is not used in calculating 
the weighted means. We define s = to be 
the center of the rotation curve, and the ex- 
tent S of each side of the rotation curve is de- 
fined by outermost s« which satisfies the S/N 
limit. We then perform a weighted mean on 
each side for all points with S/2 < Sj < S. 
Courteau (1992) obtains his best internal 
consistency with a similar algorithm. 

4.4. Weighted Mean Method: Half- 
light Radii 

We improve upon the "outer segments" al- 
gorithm by averaging the V; L (taking a weighted 
mean) over an interval in s located a fixed 
number of half-light radii from the center of 
the galaxy. Half-light radii r 1 / 2 (actually half- 



light semi-major axes) are determined from 
the /-band photometry and thus, unlike the 
"lengths" of the rotation curves above, are 
independent of galaxy distance or rotation 
curve S/N. We define the velocity of each 
arm to be the weighted mean over the range 
[ri/2, 2ri/ 2 ]- We also calculate widths us- 
ing the range [ri/2, 4ri/2]. For comparison 
we note that an exponential disk galaxy has 
r 1 / 2 equal to 1.68 scale lengths, and also has 
?"i/2 = 0.52i? op t, where i? op t is the 83%-light 
radius defined by Persic, Salucci, and Stel 
(1996, hereafter referred to as PSS96). 

4.5. Fixed-Point Method 

We measure the velocity of each arm of the 
rotation curve at a fixed number of half-light 
radii from the center. The Vi are interpo- 
lated to the appropriate point; only points 
with S/N > 3.5 are used for interpolation. 
We find that the lowest scatter is obtained 
by choosing the fixed point at 1.3ri/ 2 from 
the center, though the TFR scatter depends 
only weakly upon the exact choice as long as 
the fixed point is beyond ri/ 2 . Giovanelli et 
al. (1997b), for example, choose to measure 
W at i? op t = l-9ri/ 2 . We choose not to use 
such a large radius because for a few of our 
galaxies (and those of Matthewson, Ford, and 
Buchhorn [1992] which Giovanelli et al. use), 
this requires extrapolation past the high-S/N 
region of the rotation curve. 

4.6. Probable Min-Max Method 

The above algorithms determine some av- 
erage rotation velocity, or velocity on a ro- 
bustly determined section of the rotation curve 
It may be that the TFR is tighter when W is 
a measure of the maximum rotation speed of 
the galaxy instead of some average rotation 
speed. Unfortunately a simple maximum of 
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the Vi is not robust (Sec. 4.1). Here we define 
a fairly robust means of quantifying the ex- 
tremes of the rotation curve. We define f max 
as that which satisfies 



>Vi) = 0.1; 



(1) 



where the probabilities P are given by the 
standard error function if we assume that V; L 
has a gaussian distribution about its mea- 
sured value, with dispersion o~ Vi . In effect we 
ask the following question: At what f max is it 
likely (at the 90% level) that some part of the 
rotation curve exceeds f max ? A single high 
point on the rotation curve with a large mea- 
surement uncertainty does not unduly influ- 
ence this estimator, and it is fairly insensitive 
to binning as well. We define an analogous 
v min on the other side of the rotation curve 
and define W as the difference i> max — i> m in- 

5. Tully-Fisher Scatter 

For each of the above algorithms for mea- 
suring W we evaluate the slope and scatter of 
the resultant TFR. As explained in Paper I, 
we fit to the TFR data a model of the form 



/tot - 5 log(^/0.0233) 
-(1- ft/2 + ifej)(z- 0.0233) = 
/ + mlog( 400 ^ s _ i ) + a f e. 



(2) 



The left-hand side is the apparent total mag- 
nitude corrected to a standard redshift of z — 
0.0233, or a recession velocity of 7000 km s -1 , 
assuming all galaxies to be in free Hubble ex- 
pansion. Recession velocities are taken as the 
mean of the two arms of the rotation curve, 
and in general agree very well with the radio 
data. The /c-correction factor kj is fixed at 
0.6. On the right-hand side, e is the galaxy 
ellipticity as determined from the /-band sur- 
face photometry, and W is the measured W 



corrected to edge-on by the sine of the incli- 
nation angle, as described in Paper I. Note no 
correction is made for turbulent motion. The 
free parameters in the fit are I®, m, and the 
inclination correction coefficient aj. 

Table 1 lists the reduced photometric data 
for the 25 galaxies (raw data may be found 
in Paper I) along with the rotation and re- 
cession velocities determined from the optical 
and radio data. In Table 2 we show the results 
of the TFR fits using widths derived by the 
above algorithms. Figure 2 shows an example 
Tully-Fisher plot produced with the Probable 
Min-Max widths. 

The TFR scatter based on optical W data 
is found to be quite robust. Nearly all meth- 
ods for determining W give an RMS TFR 
scatter of ~ 0.14 mag. Even the simple Min- 
Max method works well if we consider only 
points with S/N > 3.5 (see §4.1). The only al- 
gorithm noticeably less accurate than the oth- 
ers is the Entire Side Weighted Mean, which 
leads us to believe that one should avoid the 
central portions of the rotation curve in de- 
riving W. 

5.1. Intercomparison of Widths 

Although it is not our goal to produce 
optical widths which best match the 21-cm 
widths, we do note that there are system- 
atic differences among the different optical 
W estimators, and between optical and 21- 
cm estimates of W. As expected, those al- 
gorithms which attempt to measure the ex- 
tremes of the rotation curve tend to give 
larger W values than those which attempt to 
measure the mean over some segment of the 
rotation curve. This can be seen by compar- 
ing columns (8) and (9) in Table 1. In galaxies 
which are more face-on or have smaller rota- 
tion speeds, we expect the turbulent motions 
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Fig. 2. — Tully-Fisher relation for the 25 galaxies using widths derived from Ha rotation curves 
with the "Probable Min-Max" method described in the text. Widths have been corrected to 
edge-on as described in Paper I, and /-band total magnitudes have been corrected to a common 
distance of 7000 kms" 1 , and corrected for inclination using the value of a/ shown in Table 2. The 
RMS deviation from the best-fit dashed line is 0.14 mag. 



to be a larger fraction of the observed veloci- 
ties. If different width measurement methods 
have varying susceptibility to turbulent mo- 
tions, this might be manifested as changes in 
the slope m and/or inclination correction aj 
in the derived TFRs. Note, though, that the 
values of m and a/ in Table 2 are consistent 
among various algorithms, and are consistent 
with the values obtained using 21-cm widths. 

Figure 3 plots the ratio of Probable Min- 
Max width to 21-cm width for all 25 galaxies 
function of inclination-corrected width 
Wq. No trend with Wo is present; the optical 
widths are, on average, 0.96 times the 21-cm 
widths, with no detectable dependence upon 
Wq or upon inclination. For this sample of 
relatively large, edge-on galaxies, the effect of 
turbulent velocities upon the 21-cm profiles 
(or perhaps the effect of extinction on the Ha 



rotation curves) is small, or at least constant 
across the full sample. 

This 25-galaxy sample also gives us little 
guidance on how to best treat rising rotation 
curves. In Figure 1 we see that several of the 
galaxies have rotation curves which are still 
rising at the ends, particularly C19. Should 
we take the "true" width to be the largest 
observed width (using some form of Min- 
Max method), or something larger (an ex- 
trapolation), or something smaller (one of the 
weighted mean algorithms) to get the lowest 
TFR scatter? Which of the algorithms pro- 
duces the best agreement with 21-cm widths? 
The near-equality of TFR scatter for all our 
algorithms means that we cannot offer advice 
on this question. The location of C19 near 
the top of Figure 3 suggests that the probable 
Min-Max W estimate of the rising rotation 
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Fig. 3. — Ratio of the velocity widths derived 
from Ha rotation curves (using the "Prob- 
able Min-Max" method) to the 21-cm ve- 
locity widths from Paper I, plotted vs the 
inclination-corrected optical width. Optical 
widths are on average a fixed fraction (96%) 
of the radio widths. There is no detectable 
correlation between this ratio and the magni- 
tude of the width, nor with disk inclination 
(latter not illustrated here). 

curve overestimates the rotation speed rela- 
tive to 21-cm, so an extrapolation of the ris- 
ing curve would make agreement worse. Even 
so, C19 is not an outlier from the TFR in Fig- 
ure 2. The one outlier in the optical vs 21-cm 
comparison, C9, is indeed the furthest outlier 
in the optical TFR (but is not an outlier in the 
21-cm TFR). In fact most of the extra scatter 
in the optical TFR over the 21-cm scatter is 
attributable to C9. The rotation curve of C9 
is, however, beautifully flat and square, so a 
change of algorithm would not lessen its de- 
parture from the optical TFR. Furthermore, 
C9 has no unusual features which would jus- 



tify excluding it from our sample. 

6. Rotation Curve Shapes and Third 
Parameters 

In principle the rotation curve provides 
more information than a simple width W— 
it has been suggested that some quantifica- 
tion of the shape of the rotation curve could 
be used to decrease the TFR scatter (Persic 
& Salucci 1990). Large homogeneous sam- 
ples of rotation curves show that fast-rotating 
galaxies tend to have rotation curves which 
rise more steeply in the center than slowly- 
rotating galaxies, with the fastest-rotating 
galaxies actually having a peak in the rota- 
tion curve (PSS96). 

We wish to test whether, in our limited 
sample, there is a trend of rotation curve 
shape vs W , or whether rotation curve shape 
can be used as a third parameter to reduce the 
TFR scatter. For each galaxy we compare the 
under both arms of the rotation curve to 
the area under a perfectly flat rotation curve: 

f = SUm ' 1 ' — ■ (3) 

Here s is the distance from the center of the 
galaxy, v is the recession velocity at the cen- 
ter, and ±si im are the outer bounds of the 
integration. The rotation velocity v rot is one- 
half the width W, so the value of / depends 
upon our choice of width algorithm. If we 
choose v rot = \v(+si im ) — v(—s\ iin )\/2, a solid- 
body rotator, would have / = 0.5, while a 
perfectly flat rotation curve would have / = 
1. 

Table 1 contains the calculated value of 
shape parameter / for each galaxy. We de- 
termine v r ot using the Weighted Mean method 
over the range l-4ri/ 2 , and set s\ im = 1.3ri/2- 
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Fig. 4. — Rotation curve shape vs rota- 
tion speed for the 25 galaxies. The rotation 
curve shape, parameterized by the quantity 
/ defined in Equation (2), shows no signifi- 
cant correlation with inclination-corrected ro- 
tation width Wq on a galaxy-by-galaxy basis. 
A galaxy with solid-body rotation curve in- 
side 1.3ru 2 would have / = 0.5, while a step- 
function rotation curve would have / = 1. 

For comparison with the study of PSS96, we 
note that they define R opt to be the radius 
enclosing 83% of the light of a galaxy. For 
an exponential disk with scale length R sca \ e , 
R opt = 3.2i? scale , while r 1/2 = 1.68i? scale , so 
-Ropt = l-9ri/2- Thus, / should be closely re- 
lated to the "inner slope" defined in PSS96, or 
to the rotation curve slopes defined by other 
investigators (cf. Vogt 1994 and references 
therein) . 

In Figure 4 we plot / vs Wq for the 25 
galaxies. There is no significant correlation, 
indicating that, at least in this sample, any 
trend of rotation curve shape vs amplitude is 
too weak to be detectable in individual rota- 



tion curves. This is, however, still consistent 
with PSS96, because: (1) we are looking at 
only a 2 mag range in Mj, and the PSS96 
"universal" rotation curves do not change 
drastically over such a range, and (2) even 
in the PSS96 data, the shape vs W corre- 
lation is clearest only after co-adding many 
rotation curves, and would not likely be de- 
tectable for 25 galaxies. Note that galaxies 
Cll and C17 have small / values although 
substantial parts of their rotation curves are 
flat. This is because the flat parts lie outside 
our integration range of ±1.3ri/2| a broader 
integration range would extend into regions 
of poor S/N for some of our rotation curves 
and thus require undesirable extrapolation. 

We also find that adding / parameter 
to the TFR fit in Equation || does not reduce 
the resultant scatter. Thus knowledge of the 
steepness of the rotation curve does not, in 
this sample, seem to improve our ability to 
predict absolute I magnitude. 

7. Ramifications 

We find that the RMS scatter in the TFR 
using widths from long-slit Ha rotation curves 
is ~ 0.14 mag in this sample. The TFR slope 
and scatter are robust under choice of algo- 
rithm for extracting the width W from the 
rotation curve, as long as we ignore the cen- 
tral section of the rotation curve and do not 
try to consider single points with large uncer- 
tainties (S/N < 3.5 or a v > 12kms _1 ). Using 
F-test statistics we find that the optical scat- 
ter is greater than the 0.10 mag scatter in the 
21-cm TFR for these galaxies at 95% confi- 
dence level. If we ascribe all of the 0.14 mag 
scatter in the optical-PU-based TFR to errors 
in the W's, the TFR slope of m ~ —5 im- 
plies an RMS error of 6% in the typical W 
measurement. Repeat measurements of a few 
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galaxy rotation curves have been made, some 
with slit position angle intentionally placed 
10°-20° away from the major axis, and the 
resultant W values are found to vary by 5% 
or less; we believe that measurement error in 
W is not the dominant source of scatter in the 
optical TFR. A fit to the inverse Tully-Fisher 
relation gives nearly identical results because 
the scatter is quite small. 

"Errors" in the W values can take two 
forms. First, there are measurement errors, 
which would cause the measured value of W 
to differ each time we observed the galaxy. 
Such errors might include those due to pho- 
ton statistics or poorly chosen slit position 
angles. Second, there could be what we will 
term "astrophysical" errors — reasons why the 
Ha widths do not properly measure whatever 
physical quantity is the underlying basis of 
the TFR. Examples of such cases might be 
velocities of Hn regions which depart from 
the actual rotation velocity of the galaxy, or a 
poor choice of algorithm for measuring W so 
that we do not measure the relevant rotation 
speeds properly. Of course since we do not 
actually know what the underlying physical 
basis of the TFR is in any detail, we cannot 
address these possible errors. 

The scatter in our 21-cm TFR for these 
25 galaxies is extraordinarily low; most other 
TFR surveys find RMS variations of 0.3- 
0.4 mag (Paper I). This could be low by a 
statistical fluke: our draw from the pool of ei- 
ther measurement errors or astrophysical er- 
rors (i.e. choice of galaxies) has been lucky. 
In either case, the low scatter found with 
Ha widths means that any additional errors 
caused by the use of optical widths is small 
(6% or less in W). Thus for practical pur- 
poses, we conclude that high-quality Ha ro- 
tation curves may be used in place of 21-cm 



data for TFR studies, with no practical loss 
in precision. This is also evident in the study 
of Willick et al. (1996), who intercompare 
various TFR surveys. They derive RMS scat- 
ters of 0.40, 0.38, and 0.47 mag for three sur- 
veys which use 21-cm widths, and 0.38 and 
0.43 mag for two surveys based primarily on 
Ha widths. 

The scatter in the 25-galaxy sample could 
be very low because of some as- yet-undiscovered 
physical commonality that is not present in 
the larger TFR surveys (as the lower slope of 
our TFR might also suggest). In this case we 
could also conclude that the measurement er- 
rors on the optical W values are small, a result 
which is probably extendable to larger sur- 
veys. If these 25 spirals differ physically from 
the overall spiral galaxy population, then our 
results on the robustness of the W estimators 
and on the insignificance of rotation curve 
shape may not be applicable to the general 
spiral population. 

We note also that the TFR slope and the 
inclination correction coefficient ai are ro- 
bust, with all optical W algorithms giving 
values consistent with the 21-cm values. The 
constancy of the TFR slope implies that Ha 
widths are in fixed proportion to the 21-cm 
widths over the range in galaxy sizes explored 
here. Note that our sample consists primar- 
ily of large spirals. The consistency of ai im- 
plies that none of the algorithms is excessively 
subject to inclination effects. This means, 
for example, that the measured Ha widths 
are not significantly changed relative to the 
"true" widths by the effects of extinction in 
nearly edge-on galaxies. Since we made no 
corrections for turbulent velocities in either 
the 21-cm or the Ha data, we also infer that 
such corrections are not too important for the 
high-inclination, high-Wo galaxies in our sam- 
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pie. 

Finally we note that we were unable to ex- 
tract any further useful information (for TFR 
purposes) from the rotation curves beyond 
the width W. A crude measure of rotation 
curve shape did not improve the TFR fit, nor 
did we find shape to be correlated with W 
(and presumably galaxy mass). While such a 
correlation may be present when large num- 
bers of galaxy rotation curves are co-added, 
or if a broader range of intrinsic galaxy sizes 
are studied, the galaxy-to-galaxy variations 
seem large enough to mask the effect for our 
individual objects despite the high S/N of our 
rotation curves. 
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Table 1 
Galaxy Data 
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Note. — Columns (1) and (2) give our internal name and the Zwicky or UGC name of each galaxy, 
respectively. Column (3) is the total /-band magnitude, corrected to a distance of 7000 km s^ 1 as per the 
left-hand side of Equation (1). Column (4) gives the ellipticity determined from the /-band surface photometry 
(see Paper I), column (5) is the position angle used for the long-slit observations, and column (6) is the semi- 
major axis of the isophote enclosing half the /-band light. The recession velocity in column (7) is determined 
from the optical data. Widths in column (8) are from the Probable Min-Max method, in column (9) using 
the Weighted Mean method over the range 1 2r 1 / 2 , and in column (10) we repeat the published 21-cm width 
(Paper I). Column (11) gives the rotation curve shape parameter defined in Equation (2). 
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Table 2 

Tully- Fisher Results for Various Algorithms 



Algorithm 


m 




RMS scatter 




[mag] 


[mag] 


[mag] 


Min-Max 


-5.41 ± 0.26 


1.46 ± 0.20 


0.140 


Weighted Mean: Entire Side 


-3.56 ± 0.22 


1.33 ± 0.25 


0.175 


Weighted Mean: Outer Half 


-5.47 ± 0.26 


1.35 ± 0.20 


0.138 


Weighted Mean: 1-2 


-4.95 ± 0.24 


1.18 ± 0.20 


0.139 


Weighted Mean: l-4r 1/2 


-5.08 ± 0.22 


1.25 ± 0.18 


0.128 


Fixed Point: 1.3ri/ 2 


-5.00 ± 0.24 


1.45 ± 0.20 


0.141 


Probable Min-Max* 


-5.40 ± 0.25 


1.47 ± 0.25 


0.135 


21-cm Widths 


-5.61 ± 0.18 


1.42 ± 0.13 


0.094 


^This algorithm is used to produce Figures 


2 and 3. 
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